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We study the phase behavior of ternary amphiphilic systems in the framework of a curvature 
model with non- vanishing spontaneous curvature. The amphiphilic monolayers can arrange in differ- 
ent ways to form micellar, hexagonal, lamellar and various bicontinuous cubic phases. For the latter 
case we consider both single structures (one monolayer) and double structures (two monolayers). 
Their interfaces are modeled by the triply periodic surfaces of constant mean curvature of the fami- 
lies G, D, P, C(P), I-WP and F-RD. The stability of the different bicontinuous cubic phases can be 
explained by the way in which their universal geometrical properties conspire with the concentration 
constraints. For vanishing saddle-splay modulus R, almost every phase considered has some region 
of stability in the Gibbs triangle. Although bicontinuous cubic phases are suppressed by sufficiently 
negative values of the saddle-splay modulus R, we find that they can exist for considerably lower 
values than obtained previously. The most stable bicontinuous cubic phases with decreasing k: < 
are the single and double gyroid structures since they combine favorable topological properties with 
extreme volume fractions. 
PACS numbers: 61.30.Cz, 64.70.Md, 68.10.-m, 82.70.-y 
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I. INTRODUCTION 



As a function of concentrations and temperature, amphiphilic systems generically form many different pliases, 
each of which corresponds to a specific geometrical arrangement of the amphiphilic interfaces [y|2|. Most ternary 
amphiphilic systems feature the disordered micellar, the hexagonal and the lamellar phase. Here the amphiphilic 
monolayers, which separate regions of water from regions of oil, form spheres, cylinders and lamellae, respectively. 
However, many of these systems also have stable cubic phases fSj-jg]. Near the disordered micellar phase one often 
finds a micellar cubic phase and near the lamellar phase a bicontinuous cubic phase. In the micellar cubic phase, 
the interfaces form spheres as they do in the disordered micellar phase, only that now they are packed in an orderly 
fashion. In the bicontinuous cubic phase, they form sheets which span the whole sample in all directions of space 
[^,0. Each of these triply periodic surfaces has a cubic Bravais lattice and divides space into two unconnected but 
intertwined labyrinths filled with water and oil, respectively. 

In this work, the polymorphism of ternary amphiphilic systems is studied in the framework of a curvature model with 
non- vanishing spontaneous curvature. We assume that the shape of the amphiphilic monolayers is determined by their 
bending rigidity even in the presence of concentration constraints. Therefore we model them by surfaces of constant 
mean curvature. This encompasses the spheres, cylinders and lamellae of the non-cubic phases; for the bicontinuous 
cubic phases, the surfaces of constant mean curvature are triply periodic. Safran and coworkers considered the same 
geometries and showed that the interplay between the bending energy and the volume constraints can explain some 
aspects of the generic phase behavior of ternary amphiphilic systems |P-p^. In particular, Wang and Safran |l^] 
considered the single structure (one monolayer) and the double structure (two monolayers) from the D-family. Their 
calculation relied on data for triply periodic surfaces of constant mean curvature, which were found numerically by 
Anderson, who investigated the geometrical properties of the D, P, C(P), I-WP and F-RD families However, 
the bicontinuous cubic phases most often identified experimentally are the various gyroid structures which correspond 
to the G-family. Data for the G-family were calculated only recently by Grofie-Brauckmann [^jwj . In this work, we 
consider as bicontinuous cubic phases single and double structures for all families for which the required geometrical 
data are known, including the G-family. Indeed, we find that for negative values of the saddle-splay modulus k, 
of all bicontinuous cubic phases considered it is mainly the gyroid structures which are stable. They can exist for 
values of R, which are twice as negative as the values for which the structures from the D-family have been found 
to be stable previously [l0|] . We will show that this finding can be explained nicely by considering the interplay 
between certain universal geometrical properties of the various bicontinuous cubic phases and the volume constraints 
within the Gibbs triangle. In fact we find that the gyroid structures are so favorable because they have exceptional 
topological properties and at the same time can accommodate extreme volume fractions of oil and water. 



II. CURVATURE MODEL AND NON-CUBIC PHASES 



The elastic properties of amphiphilic interfaces are described by the Canham-Helfrich expression for the elastic 
energy per unit area [^p^ 

felasUc =2k{H- C^f +RK , (1) 

where H = (ci +C2)/2 and K — C1C2 are mean and Gaussian curvatures, respectively, and Ci and C2 the two principal 
curvatures of the surface. The two elastic moduli are the bending rigidity k and the saddle-splay modulus R. The 
main temperature dependence is carried by the spontaneous curvature cq. For systems with water, oil and non-ionic 
surfactants CiEj, it is found experimentally that cq oc (Tf, — T) where Tf, is the balanced temperature Here 
positive curvature is defined to be curvature towards the oil regions, thus oil is the interior phase below the balanced 
temperature and water above. For example, cq ~ 1/(60 ^^'^ ~ 1/(12^ for H2O/C14/C12E5 at T = 20°C and 
T = 34°C, respectively, where / « 1.5 nm is the amphiphile length in the monolayer and Tf, = 48°C In this 

work, we assume a positive value for the spontaneous curvature, as it is typical for surfactant systems below the 
balanced temperature, so that oil is the interior phase. The tendency to bend towards the oil regions decreases with 
temperature since the headgroups' hydration decreases. For the case of a negative spontaneous curvature, water and 
oil have to be interchanged in the various structures as well as in the phase diagrams presented below. This is usually 
the case not only for surfactants above the balanced temperature, but also for lipid systems where the monolayers 
tend to bend towards the water regions due to their bulky tail regions. 
We consider the dimensionless free energy per unit volume. 
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where the integration extends over the neutral surface A of the amphiphilic monolayers, V is the overall volume, and 
r = —k/2k. For the first (bending) term in Eq. (|^), the area integration can be carried out since we only consider 
surfaces of constant mean curvature. For the second (topological) term, we employ the Gauss-Bonnet theorem, 
J dAK — 27TX, where x is the Eulcr characteristic of the surface. The curvature model is only stable for —2k < k < Q 
or < r < 1 jl^ . Experimentally there is no straightforward way to measure k, but usually a small negative value is 
assumed. 

The phase diagram at constant temperature is a function of the volume fractions pw, po and pA of water, oil and 
amphiphile, respectively, which are restricted to the Gibbs triangle by pw + Po + PA = 1- We consider amphiphiles 
of length /, which have a tail length al and a head size (1 — a)l, with < a < 1. Throughout this paper we will use 
a = 1/2, i.e. we will consider amphiphiles which have head and tail regions of similar size, as they are for C12E5. 
The hydrocarbon volume fraction is given hy v = pQ + apA |9|. In order to parametrize concentration space, it is 
convenient to use v and the ratio w — pa/Hpo + o:pa)coI]- The amphiphile volume fraction pA is taken to be Al/V. 
If the amphiphile is assumed to occupy the space between the two parallel surfaces at distances al and (1 — a)l from 
the neutral surface, this is a very accurate approximation when the amphiphile length / is small compared to the 
extension of the amphiphilic aggregates, and gives reasonable results even for the extreme case of oil-free spherical 
micelles for a ~ 1/2. 

In the following all lengths are measured in units of I. The free-energy densities of the non-cubic phases follow from 
Eq. (I) as 

(3) 
(4) 

(5) 

Note that only fs depends on r, since the other two structures have no Gaussian curvature. Since we only consider 
the free-energy contribution due to the curvature elasticity of the amphiphilic monolayers, for the non-cubic phases 
/ has a trivial u-dependence. However, due to close packing, spheres and cylinders have maximal volume fractions 
Vmax — V27r/6 ~ 0.74 and tt/2^/3 = 0.91, respectively. With w, v and r, the model has a three-dimensional parameter 
space. 

It can be seen from Eq. (3) that the Maxwell construction is not possible for this model, since is not concave. 
Therefore we use the intersections of the free-energy densities of different phases as an indication for the location of 
phase transitions. We want to remark parenthetically that Wang and Safran ||l^ considered the free energy per unit 
area (rather than per unit volume); this amounts to an overall factor of wv/cq (the dimensionless specific area) in the 
free-energy density of all phases. As long as the location of the phase transitions is estimated from the intersection of 
the free-energy curves, the two approaches are equivalent. However, for a calculation of two-phase regions, the use of 
the free energy per unit volume becomes essential — since in general coexisting phases will have different amphiphile 
concentrations . 

For r = (vanishing saddle-splay modulus) and v < 0.74 (no excluded volume effects), Eqs. (3), (4) and (5) 
imply the phase sequence L C S with decreasing w. fs begins to rise again for w < 3. This is identified 
with the emulsification failure, the coexistence of S and an excess oil phase at low amphiphile concentration JlOt . In 
this paper, we define the emulsification failure not as the minimum in fs in regard to w, but by using a Maxwell 
construction in regard to w between S and an excess oil phase, i.e. by solving dfs/dw — fs/w for w. With increasing 
r (and v < 0.74), spheres become more favorable, while cylinders are increasingly suppressed and finally disappear 
for r > 1/4. In particular, the phase boundaries S — C, S — L, C — L and the emulsification failure are obtained from 
Eqs. (3), (4) and (5) to be w = 24/(7 - 16r), w = 6/(1 - r), w = 8 and w = 3/(r - 1), respectively. 

It is important to realize that not all values of {v, w) are physically relevant. In Fig. n we show the mapping between 
the (i;,w)-plane and the Gibbs triangle for ternary mixtures. The binary limits 0-W, W-A and A-0 of the Gibbs 
triangle correspond to po + pw = li Pw + PA = ^ and pA+ Po = 1, respectively. In the {v, w)-plane, this corresponds 
to the hues w = 0, w = 1/(q;co) and w — (1/u — !)/((! — ct)co), respectively. The line w = is mapped onto the 
W-apex. Fig. |l| demonstrates (for a — 1/2) this mapping of the {v,w)-pla.ne onto the Gibbs triangle for cq = 1/6 
and Co = 1/12. As mentioned above, these values correspond to the system H2O/C14/C12E5 a,t T — 20°C and 
T — 34° C, respectively [|l^. Smaller/larger values for the spontaneous curvature (which corresponds to higher/lower 
temperatures or other components) increases/decreases the size of the relevant region in Fig. |l|a to larger/smaller 
values of w. One also can see in Fig. 0b that for a = 1/2 lines of constant v are perpendicular to the W-0 side and 
lines of constant w are straight lines through the W-apex. 



JL(W, V) = wv, 

fc{w,v) = wv 
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fs{w,v,r) = wv 
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III. BICONTINUOUS CUBIC PHASES 



A. Properties of constant-mean-curvature surfaces 

We model the amphiphilic monolayers in the bicontinuous cubic phases by triply periodic surfaces of constant 
mean curvature H — [ci + C2)/2. For the special case H = 0, this leads to triply periodic minimal surfaces with 
ci — — C2 and K = C1C2 < 0. Due to the Gauss-Bonnet theorem, J dAK — 27rx, these surfaces have a negative Euler 
characteristic x P^r unit cell. Triply periodic minimal surfaces are commonly used to model all kinds of extended 
sheet-like structures in condensed matter systems, in particular the mid-surfaces of the lipid bilayers in inverse cubic 
phases which are very prominent in lipid- water mixtures Before 1970 only three cubic triply periodic minimal 

surfaces have been known (D, P and C(P)) [|l|,|0). Then Schoen described five more (G, F-RD, I-WP, 0,C-TO 
and C(D)) |2l]. Today some more examples are known |^,^, but none of them seems to be of physical relevance. 
Karcher proved in 1989 not only the existence of the triply periodic minimal surfaces described by Schoen, but also 
that the simpler of them can be deformed into triply periodic surfaces of constant (non-zero) mean curvature [ p4[ . 
As for any triply periodic surface, space is divided into two percolating labyrinths. A shift of H to positive/negative 
values shrinks/expands one labyrinth, while it expands/shrinks the other. Thus two branches are generated, which 
both end in cubic arrangements of (infinitcsimally connected and possibly self-intersecting) spheres. Since the Euler 
characteristic is a topological quantity connected to the genus g of the surface by x = 2(1 — g), it does not change 
within a family. Anderson studied the cubic families D, P, C(P), I-WP and F-RD and calculated for a conventional 
unit cell of unit lattice constant both the volume fraction v of one labyrinth as a function of scaled mean curvature 
H* (the volume fraction of the other labyrinth follows as \ — v) and the scaled surface area A* as a function of v 
Recently v{H*) and A*{v) were calculated also for the G-family by GroBe-Brauckmann p3| , p^ . 

In this work we consider the families G, D, P, C(P), I-WP and F-RD, for which these geometrical data are available. 
Their minimal surface member is shown in Fig. g for D, C(P), I-WP and F-RD in one conventional unit cell. The 
numerically calculated data points for v{H*) and A*{v) per conventional unit cell arc taken from Rcfs. |pT| , |l3| ; 
rearrangement and interpolation with cubic splines provides smooth functions v{H*), H*{v) and A*(v). The data are 
used only up to the extremal values of v, beyond which the curves bend backwards. Beyond these points, the surfaces 
resemble ensembles of spheroidal regions connected by nearly unduloidal necks, which we do not consider to be of 
physical relevance [^. The functions H*{v) and A*{v) as used in this work are plotted in Fig. ^ for the six families 
considered. In fact, only one of the two curves for each family carries independent information since the other one 
can be constructed by using dA* — 2H*dv ||l^. Each family only exists over a certain range of volume fractions; the 
extreme cases are the G-family with v G [0.056,0.944] and the C(P)-family with v S [0.481,0.519]. If the minimal 
surface member of a family divides space into two congruent labyrinths, they have the same volume fraction vq = 1/2, 
the two branches are symmetric to each other and the minimal surface as well as the family itself is called balanced. 
For the structures considered here, this is the case for G, D, P and C(P). For these structures, the curves H*{v) and 
A*{v) are symmetrical with respect to v = 1/2. For the non-balanced families I-WP and F-RD, the two labyrinths 
are of different topology and vq 7^ 1/2. In Tab. |, we collect data connected to the minimal surface members of each 
family, which we will use later in our discussion. 

In the vicinity of the minimal surface, the volume fraction w as a function of scaled mean curvature H* can be 
approximated as v{H*) ^ vq ~ cH* + 0{H*'^). With dA* ^ 2H*dv it follows that 

H*iv)^J-!^^ + 0{{v-vor) (6) 

A*{v)=Ao-^^^^^^+0{iv-vor) 
c 

near the minimal surface. The values for Aq, vq, c are given in Tab. |. Since the surfaces of constant mean curvature 
can be expected to have similar geometrical properties as parallel surfaces in the vicinity of a minimal surface, the 
magnitude of c can be estimated as follows. If t denotes the perpendicular distance from the minimal to its parallel 
surface, to lowest order in t the volume fraction v and the mean curvature H* , averaged over the surface in the whole 
unit cell, are given hy v = vq + Agt and H* — 27rxt/^o, respectively [ p^ , p7[ |. Thus H* — 2ttx{v ~ i'o)/^o^ and 
c = — Ao^/27rx for the parallel surface case. The corresponding numbers are given in Tab. ||; except for C(P), the 
overall agreement with the numerical data for c is remarkably good. 
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B. Single and double structures 



The simplest case of a cubic bicontinuous phase in a ternary amphiphihc system is a single structure where the 
amphiphihc monolayers form one triply periodic surface. Then one labyrinth is filled with oil and the other with 
water. For balanced families, filling either of the two labyrinths with oil gives the same single structure; for the 
non-balanced families I-WP and F-RD, this yields different single structures which we denote by I, WP, F and RD, 
respectively Here the symbols I and F correspond to the curves plotted for I-WP and F-RD in Fig. ||. Since 

mean curvature is defined as curvature towards the oil regions, WP and RD follow by the replacements v 1 — v 
and H ^ —H from the data of I and F plotted in Fig. respectively. Thus, altogether we consider 8 different 
single structures, which exist for the volume intervals [0.056,0.944] for G, [0.131,0.869] for D, [0.249,0.751] for P, 
[0.481,0.519] for C(P), [0.357,0.857] for I, [0.143,0.643] for WP, [0.439,0.625] for F and [0.375,0.561] for RD. 

To each single structure corresponds a double structure where the amphiphihc monolayers form two triply periodic 
surfaces arranged roughly parallel and on either side of the minimal surface of the corresponding single structure. Both 
surfaces divide space into two labyrinths which are topologically equivalent to those of the initial structure. However, 
since two amphiphihc monolayers are present, now these labyrinths are filled with the same component and separated 
by a bilayer which is filled with the other component. Single and double structures are also known as monolayer and 
bilayer structures [|2^,^ . A double structure can be either of type I (water-filled bilayer) or of type II (oil-filled bilayer) 
[ p8[ . There will be no problem below to distinguish the symbol I for double structure of type I from the symbol I for 
one of the two simple structures of the I-WP family. In Tab. |l|, we summarize the classification of bicontinuous cubic 
structures. In the following, single structures, type I double structures and type II double structures are abbreviated 
as Q, Q/ and Q//, respectively. Q// structures are also known as inverse bicontinuous cubic phases. In Fig. ^ we 
show single and double gyroid structures in one conventional unit cell. For many surfactant-water and lipid-water 
systems, G// and/or D// are well established; moreover there are reports on P// and G/ structures in these binary 
amphiphihc systems [||-||]. Q structures in binary systems have been discussed only theoretically so far [29|-|3l|. For 
many systems with water, oil and surfactant, stable bicontinuous cubic phases have been reported, however often 
without identification of their space group (see e.g. Ref. for the system H2O /Ciq/Ci2E^). The best established 
identification is a number of Q// structures for the system DDAB-water-styrene There are a few reports on Q 
structures |l^,Q, and none on Q/ structures. However, several speculations on single structures can be found in 
the literature (e.g. in Refs. |27||3^ ]), and our recent theoretical work on ternary systems with vanishing spontaneous 
curvature suggests that they should have some physical relevance [^^. To our knowledge, hardly nothing is known 
on bicontinuous cubic phases in ternary systems with water, oil and lipid. 

In order to construct a double structure from a given family of triply periodic surfaces of constant mean curvature, 
we take two surfaces corresponding to H and —H. Thus there is one Q/ and one Q// structure for each of the 
6 families. For Q/ structures, the minimal surface case corresponds to w = 1. As both surfaces accumulate mean 
curvature in their respective branches, the volume fraction decreases until the first labyrinth reaches its minimal size 
(which is the minimal volume fraction of the corresponding single structure). The volume intervals covered by the Q/ 
structures turn out to be [0.112, 1.0] for G/, [0.262, 1.0] for D/, [0.498, 1.0] for P/, [0.962, 1.0] for C(P)/, [0.624, 1.0] for 
I-WP/ and [0.818, 1.0] for F-RD/. The Q// structures follow from the corresponding Q/ structures by an interchange 
of oil and water. Thus here the minimal surface case corresponds to w = 0, and the volume fractions covered are 
complementary to the ones given for the Q/ structures. 



C. Curvature and topology index 



For bicontinuous cubic structures, the free-energy density depends on the hydrocarbon volume fraction u in a non- 
trivial way, since only one or two amphiphihc aggregates are present. The lattice constant is denoted with a. For a 
given value of v, the mean curvature H{v,a) — H*{v)/a and the surface area A{v,a) = A*{v)a?' within a unit cell 
are determined by the curves plotted in Fig. ^j. The amphiphile concentration pA = A{v, a)/a^ — A*{v)/a fixes a, so 
that a = A*{v)/pA — A* (v)/ (wvcq). Using the Gauss-Bonnet theorem for the Gaussian curvature term, wc find from 
Eq. (^ for the free-energy density of a bicontinuous cubic structure 

r . .A (H*{v) X 27Txr 
jBc{w,v,r) = — - 
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wv [A[v) wv — 1\ + r p^^p ^ > 
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with 



H*{v) 




,3\ 1/2 



A(^;) 



A* {v) 



(8) 



Here x denote the Euler characteristic of the interfaces within one conventional unit cell, as given in Tab. ^ The 
curvature index A.{v) and the topology index T{v) are two universal geometrical quantities which characterize a surface 
in three-dimensional space. Their significance can be understood rather easily from the observation that both the 
mean and Gaussian curvatures can be made dimensionless by multiplying them with appropriate powers of V/A, 
which is the only relevant length scale; this implies H{V/A) = A and |iir|(F/A)^ = 2Tr\x\VyA^ = It is not 

difficult to see that both quantities are not only invariant under scale transformations, but also under a change of unit 
cell. They also occur in integral geometry, a mathematical theory concerned with invariant geometrical measures [ p7[ . 
For convex bodies in three-dimensional space there exist two independent isoperimetric inequalities and therefore two 
independent isoperimetric ratios, which usually are chosen to be A and l/F^. The equivalence between these quantities 
and the isoperimetric ratios holds since for surfaces of constant mean curvature, the integral mean curvature, which is 
one of the Minkowski functionals of integral geometry, can be replaced by HA. The curvature index A describes how 
strongly the structures is curved and the topology index F describes its porosity (the larger its value, the less holes 
the structure has). For minimal surfaces, A vanishes and F remains the only relevant quantity. Its significance for 
amphiphilic systems has been pointed out by Hyde pTf , and its variation as a function of crystallographic determinants 
has been studied recently by Fodgen and Hyde js^. For multicontinuous structures build from n sheets, the curvature 
index is smaller and the topology index is larger by a factor of n compared to the corresponding single structure. 
In particular, for double structures the curvature index is half and the topology index is twice as large as for single 
structures. Thus, double structures are less porous than single structures since they have disconnected surfaces. We 
can infer from Tab. ^ that the minimal gyroid is the least porous of the single structures with F = 0.7667. We want to 
remark parenthetically that — surprisingly — bicontinuous random surfaces have similar values of F as single cubic 
structures. For example, it can be shown exactly that the isosurfaces of Gaussian random fields with (H) = have 
F = Vs/n = 0.9003 ||,|o|. Since this value is only slightly larger than that of the single gyroid phase, it can be 
concluded that the random sponge phase on average features only few disconnected or multiple sheets. However, 
when the random sponge's interfaces are made to acquire curvature, the topology index grows strongly with A since 
disconnected parts proliferate pO[| . The opposite is true for the families of surfaces of constant mean curvature 
considered here; if the mean curvature is increased from zero, the topology remains the same, but the area content 
decreases. Therefore the topology index decreases, since for given topology it measures the specific surface area. 

In order to evaluate the free-energy density of the bicontinuous cubic phases from Eq. (7) as a function of w, v 
and r, curvature and topology index as defined in Eq. (||) have to be calculated as a function of v. For the 8 Q 
structures, they follow in a straightforward way from the data plotted in Fig. |[ In order to derive them for the 6 
Qi structures, we first consider the non-balanced case, thus the two labyrinths 1 and 2 have different topologies (e.g. 
for I-WP, we have 1 = 1 and 2 = WP). We first construct vi{H*) = vi{H*) + V2{H*) and then invert it to obtain 
Hj{v), the mean curvature of both surfaces as a function of the overall volume fraction (which is distributed onto 
both labyrinths). The surface area then follows as A*j{v) = Al{vi{Hj {v))) + A2{v2{Hj {v})). The Euler characteristic 
is X/ = Xi + X2 (with xi = X2 given in Tab. ||). With Hj{v), A*j{v) and xi calculated, we then can evaluate A and 
F in Eq. (||) for the non-balanced Q/ structures. For the non-balanced Q// structures, we have to evaluate them 
using Hjj{v) = —H'j{l — v), A*jj{v) = Aj{l — v) and xii ~ Xi- For the balanced families, one has vi = V2 = w/2, 
H*j{v) = H*{v/2), A*j{v) = 2A*{v/2), H*ji{v) = -H* {{I ~ v) /2.) A*jj{v) = 2A*((1 - v)/2). This amounts to using 
Ai{v) = A{v/2)/2, F/(u) = 2F(u/2), Aji{v) = -A((l - v)/2)/2 and Tii{v) = 2F((1 - v)/2) in Eq. §). Note the 
minus sign for the curvature of the Q// structures since the monolayers change their orientation compared to the Q/ 
structures. In Fig. || we plot A and F as a function of v for all single and double structures considered. Recall that 
the minimal surface case corresponds to v ~ vq, 1 and for Q, Q/ and Q// structures, respectively; in these cases, 
the curvature index disappears and the topology index acquires its maximum value. In Tab. |, we give the values of 
Fq = F(w = uq) for all families considered. This implies the hierarchy G, D, I-WP, P, C(P), F-RD with decreasing 
Fq. From Fig. ^ we also see that the single structures I and WP in fact can become better than single D for certain 
values of v. 

We see from Fig. ^ that single structures for v > vq and all Q// structures have negative values for the curvature 
index A. From Eq. (7) it therefore follows that these structures have fsc > wv = Jl, compare Eq. (3); thus they are 
always less stable than the lamellar phase. In other words, since we consider the case that the monolayers prefer to 
bend towards the oil regions (positive spontaneous curvature cq), in the framework of the curvature model no phase 
can be stable which curves towards the water regions. One should note, however, that single structures for v > vq 
or Q/7 structures could occur even for positive cq if they were stabilized by other contributions to the free energy 
which are neglected in our treatment. For the case of negative cq, the situation is reversed: now the monolayers prefer 
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to bend towards the water regions, and from the double structures the (inverse) Qu and not the Q/ structures are 
stable. In general, for ternary systems one expects Q in the middle of the Gibbs triangle (possibly in the vicinity 
of a microemulsion phase for surfactant systems), Q/ near the binary side amphiphile-oil and Q// near the binary 
side amphiphile-water. For surfactant systems, which usually have cq > at room temperature, the curvature model 
predicts the stability of Q and Q/ structures; for lipid systems, which usually have cq < at room temperature, Q and 
(inverse) Q// structures are predicted. This general prediction conforms with the predominance of inverse phases for 
lipid- water mixtures. Note that the surfactant systems DDAB-water-styrene mentioned above is an exception to 
our distinction between surfactant and lipid systems since the two tails result in cq < and therefore lead to inverse 
phases like in the lipid case. 



Altogether we consider 17 different phases, 3 non-cubic, 8 single and 6 double structures of type I. Moreover there 
always exists an emulsification failure at low amphiphile-to-oil ratios. We first discuss the case of vanishing saddle- 
splay modulus, r — 0. For the bicontinuous cubic phases we see from Eq. (7) that the optimal value / = for the 
free-energy density is achieved for w{v) = l/[vA{v)]. For the spherical and cylindrical phases we find w{v) = 3 and 
'w{v) = 4 from Eqs. (5) and (4). For each phase the corresponding line w{v) lies in the middle of its region of stability; 
these lines are plotted in Fig. ^ instead of the full phase diagram for r = 0. In the following we will denote them as 
lines of vanishing frustration since they mark the specific parameter values for which a given phase can satisfy both 
bending and concentration constraints simultaneously. In the {v,w)-plane, each of the structures considered has such 
a line, thus each of them has some region of stability where its particular geometry serves best to accommodate the 
volume fractions of the different components. In fact there are even several values of {w, v) where the free-energy 
density of different structures is degenerate. In Fig. ^ we map the lines of vanishing frustration from Fig. |^a onto 
the Gibbs triangle for cq = 1/6. For these values, the Q structures run towards the W-A side and the Q/ structures 
towards the A-0 side. Although the lines of vanishing frustration for C(P), C(P)/ and F-RD/ are not mapped onto 
the Gibbs triangle, the phase behavior here remains highly degenerate. 

Fig. ^ demonstrates that the different structural types considered occupy different regions of the phase diagram 
in a very characteristic fashion: for large w, Q and Qj structures are stable for w < 1/2 and f ^ 1, respectively. 
With decreasing w, the regions of stability curve to the left. In order to understand the sequence of phases within 
the band- like region of each structural type for large w, it is useful to expand the free-energy densities of the various 
bicontinuous structures about the minimal surface members by using Eq. (6). For the Q and Q/ structures, this 
corresponds to an expansion about the volume fractions vq and 1, respectively. Again we only consider the lines of 
vanishing frustration for which / = and w{v) = l/[vA{v)]. The curvature index A{v) defined in Eq. (^ can be 
approximated by 



In the same order of approximation, the lines of vanishing frustration then follow as w — —cAq/vo{v — vq) and 
w = —AcAo/{v — 1) for Q and Q/ structures, respectively. Thus the minimal surface case corresponds to the stable 
solutions for w ^ 1 at w = a-nd w = 1, respectively, where the curvature indexes of the corresponding structures 
disappear (compare Fig. For the balanced single structures and the double structures the hierarchy of the different 
phases within the band-like region occupied by a certain structural type is thus determined by the values of cAq, 
which are given in Tab. | for the 6 families considered. The only exception are the 4 non-balanced single structures, 
which cannot be compared in this way, since the value for vq is different for each of them. Using the approximation 
c w — Ao^/27rx derived above, we find cAq « Fq^, so that the sequence is approximately determined by the topology 
index of the minimal-surface member of each family. In fact the structures in Tab. || are ordered with decreasing 
Fg. In particular, for a given structural type and values of the hydrocarbon volume v outside the (rather restricted) 
w-intervals of existence of I-WP, F-RD and C(P), we expect from Eq. (9) to find the sequence G - D - P as a function 
of either v or w. 



IV. PHASE BEHAVIOR 



A. Phase diagrams for r — 




(9) 
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B. Phase diagrams for r > 



When the saddle-splay modulus becomes negative, so that r > 0, the free-energy densities of lamellae and cylinders 
do not change since these structures have vanishing Gaussian curvature. However, spheres have positive Gaussian 
curvature and therefore their free-energy density decreases (compare Eq. (5)). Since bicontinuous cubic phases have 
negative (integral) Gaussian curvature, their free-energy density increases (compare Eq. (7)) and they will be sup- 
pressed towards large w by the lamellar phase and towards small w by the cylindrical phase; these phases in turn 
will for sufficiently large r be suppressed by the spherical phase. Since ternary amphiphilic systems presumably have 
negative but small values of the saddle-splay modulus, the important questions here are up to which value of r the 
bicontinuous cubic phases remain stable, and which of the 14 different structures considered performs best. It follows 
from Eq. (7) that the relevant quantity is the topology index: the larger its value for a certain bicontinuous cubic 
structure which is stable for r — 0, the longer this structure stays stable with increasing r. Fig. || shows that for any 
value of V, the double structures have larger geometry indices than single structures. Within each of the two relevant 
structural classes, it is the gyroid structure which has the largest value of the topology index (compare Tab. ^. We 
therefore conclude that the double gyroid G/ should dominate phase behavior for r > for topological reasons. 
However, there are three restrictions to this general conclusion. First, G/ can only realize v S [0.112,1.0]. Second, 
before G/ can dominate all other bicontinuous cubic phases with increasing r, it might be already dominated itself by 
the lamellar phase L and the cylindrical phase C. And third, since the topological term is weighted by a factor {wv)^ 
in Eq. (7), the double gyroid cannot perform so well for small v as it can for larger v. 

Our numerical results nicely corroborate this analysis and shows the exact outcome of the balance between the 
different principles mentioned. In Fig. |^ we show phase diagrams as a function of w and r for v = 0.1, 0.2, 0.4 and 0.6. 
Since w, v and r define the parameter space, these figures show nearly the complete phase behavior predicted by the 
model. Only the subsequent mapping onto the Gibbs triangle is affected by the chosen values for a and cq. In Fig. 
we draw the lines of constant w, which marks the upper edge of the part of the (r, ?z;)-plane which is mapped onto the 
Gibbs triangle for a = 1/2 and cq — 1/6. The degeneracy of the bicontinuous cubic phases for r — discussed above 
disappears quickly with increasing r. They eventually all disappear because the lamellar phase L and the cylindrical 
phase C become more stable. At r = 0.25, C itself is suppressed by the spherical phase S. The double gyroid G/ 
remains stable for larger values of r than all other bicontinuous cubic phases. In fact for v = 0.2 (Fig. |^b) it is stable 
up to r = 0.2. This result stands in marked contrast to the results of Ref. ||l3|, which predicted that the most stable 
phase should be D/ with a stability limit of r = 0.1. For v = 0.1 (Fig. |^Jthe double gyroid cannot exist and the 
single gyroid G is the only stable bicontinuous cubic phase. For v = 0.6 (Fig. 0d) the single structures are not stable 
since their curvature index is negative for v > vq (compare Fig. ^). If more than one structure is stable within one 
structural class, we see the sequence G - D - P which was shown above to be determined by the topology index as 
well. 

In Fig. H we show the Gibbs triangles for r = 0.01 and cq = 1/6, r = 1/15 = 0.07 and cq = 1/6, as well as r = 1/15 
and Co = 1/12. In the first case of a very small value of the saddle-splay modulus, one still sees the degeneracy of the 
case r = 0. For the more negative value, only the single gyroid G and the double gyroid G/ are stable. Comparing 
with Fig. H, we see that the G-phase is stable near its extremal volume fraction of v = 0.056, while the G/-phase is 
stable for a large range of volume fractions, which covers the region of the Gibbs triangle where its line of vanishing 
frustration is located for r = 0. The effect of decreasing spontaneous curvature from cq = 1/6 in Fig. ||b to cq = 1/12 
in Fig. ^ is to extend the region of stability of the lamellar and bicontinuous phases away from the W-A side, towards 
the center of the Gibbs triangle. 

V. DISCUSSION AND CONCLUSIONS 

In order to calculate phase behavior of ternary amphiphilic systems, we investigated a simple curvature model 
with non-zero spontaneous curvature. In particular, we focused on the bicontinuous cubic phases whose interfaces 
were modeled by triply periodic surfaces of constant mean curvature. We showed that for this class of surfaces, the 
free-energy density of the bicontinuous cubic phases can be written in a very general form which emphasizes the 
universal geometrical character of the problem studied. It consists of two terms: the bending term depends on the 
curvature index A and the topological term depends on the topology index F. The relative strength of the two terms 
is determined by r = —k/2k. The relevance of the two quantities A and F not only depends on r, but also on the way 
by which they are weighted in the free-energy density by the two variables w and v which parametrize concentration 
space. Several properties have to conspire for a specific phase to be stable at a certain point of the phase diagram. 
First, its geometrical properties have to allow to accommodate the given concentrations (this imposes some constraints 
on the allowed values for the hydrocarbon volume fraction v). Second, the resulting mean curvature must be close to 
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the given spontaneous curvature in order to keep the bending term small. Since there are two independent degrees 
of freedom in concentration space, one of which is sufficient to adjust the mean curvature to its optimal value, the 
regions of stability are at least one-dimensional in the {v, it;)-plane for r — 0. The exact location of these lines of 
vanishing frustration is determined by the curvature index A{v). Phases with negative curvature index A are less 
stable than the lamellar phase. For positive spontaneous curvature (monolayer bending towards oil regions), this rules 
out the single structures for v > vq and the double structures of type II. Third, for r > the structure has to be 
favored by the topological term also. This requires large values for the topology index r{v) and basically favors the 
double gyroid G/. Fourth, the concentrations in the regions of stability have to be physically relevant, i.e. they have 
to correspond to those of the Gibbs triangle. The mapping of the phase diagram as a function of v and w onto the 
Gibbs triangle depends on specific values for the amphiphile chain length a and the spontaneous curvature cq and 
cuts off some of the phase behavior in the (w, u;)-plane. 

The phase behavior for r = is highly degenerate. This degeneracy has been discussed already in Ref. for 
the Canham-Helfrich Hamiltonian without spontaneous curvature and concentration constraints. In order to resolve 
the question of the relative stability of the different bicontinuous cubic phases, one has to consider further physical 
effects like topological contributions (r > 0), van der Waals, electrostatic or steric interactions, higher order curvature 
terms or packing energies for the hydrocarbon chains Our analysis shows that although every phase 

considered has some region of stability in phase space, there are certain general principles which allow to understand 
the complicated structure of the resulting phase diagram. The locations of the regions of stability are determined by 
the curvature index A; we have derived a simple approximation for A, which is valid for large w and explains why Q 
and structures are stable for w < 1/2 and v < 1, respectively. Moreover, it turns out that the sequence of phases 
within the band-like regions of a certain structural type is determined by the value of cAq, which in turn can be well 
approximated by Tq^. Thus the relative location of the different phases of one type is determined by the topology 
index of the minimal-surface member of that family. 

When the saddle-splay modulus becomes negative, more and more of the cubic bicontinuous phases disappear 
until eventually all of them are suppressed by the non-cubic phases. However, we found that the bicontinuous cubic 
phases remain stable for considerably higher values of r than found previously |l0| . Since the structure performs 
best which has both a high topology index and can accommodate large ranges of hydrocarbon volume fraction v, the 
double gyroid G/ becomes the most stable bicontinuous cubic phase for increasing r. Note that G/ is the only double 
structure of type I where the oil-filled labyrinths consist of channels, which are connected by junctions of three-fold 
coordination exclusively. Its outstanding stability with respect to r can therefore be explained by the fact that it is 
the most cylinder-like of the bicontinuous cubic phases. The same geometrical property in fact can stabilize also the 
entropy-dominated microemulsion since low-coordinated vertices provide a lot of configurational entropy [p| . Since 
with increasing r the double gyroid G/ is eventually suppressed by the lamellar and the cylindrical phase before it 
can suppress all other bicontinuous cubic phases, the single gyroid G (with has the most favorable topology index 
within its structural type) has a considerable region of stability as well. 

In summary, we have demonstrated that the complicated phase behavior of bicontinuous cubic phases in ternary 
systems can be understood in terms of the interplay between their universal geometrical properties and the concen- 
tration constraints of a ternary system. The main result is that the gyroid structures are favored since they have the 
largest values for the topology index F, that is the smallest porosity among the bicontinuous cubic phases. It should 
be pointed out, however, that this work rests on the basic assumption that all structures form surfaces of constant 
mean curvature. Such surfaces arise as minima of J dA under a volume constraint, or as minima of J dA{H — cq)^ 
without volume constraint, but they are not solutions to J dA{H — cq)^ under volume constraint — except for the 
special case of vanishing frustration, where the surface with H = cq just satisfies the volume constraint. It is easy to 
understand that, in general, surfaces of constant mean curvature do not minimize the curvature energy with volume 
constraint, since the energy can be lowered by keeping the mean curvature on the largest part of the surface very 
close to Co, and by concentrating deviations from H = cq — which are enforced by the volume constraint — to a 
small part of the surface (which therefore makes a small contribution to the curvature integral). However, since the 
regions of stability of most bicontinuous phases do not extend very far from the lines of vanishing frustration in the 
phase diagram, the deviations of the exact solutions from surfaces of constant mean curvature can be expected to be 
small for the physically relevant regions. 

In this work, we focused on the case of finite and positive spontaneous curvature. Its value can be controled in 
amphiphilic systems by changing temperature and disappears at the balanced temperature. Above the balanced 
temperature, the spontaneous curvature is negative and oil and water have to be interchanged in the structures and 
phase diagrams discussed. Then the double structures of type I are replaced by the double structures of type II (the 
inverse phases) and the phase boundaries in the Gibbs triangle run predominantly towards the oil apex. For surfactant 
systems, the balanced temperature is usually well above room temperature [|7|. Although the curvature model 
presented here predicts inverse bicontinuous phases above the balanced temperature, they are probably destroyed by 
thermal fluctuations in this case. For lipid systems, the balanced temperature is usually similar to the main transition 
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temperature, and inverse phases are predicted even at room temperature. In accordance with experiments on both 
kinds of systems, the model presented in this paper predicts that the lamellar phase dominates at the balanced 
temperature. However, since it assumes finite spontaneous curvature and neglects thermal fluctuations, it cannot 
describe the bicontinuous microemulsion phase which for surfactant systems often coexists with the lamellar phase 
around the balanced temperature. 

We did not consider the effect of thermal fluctuations or long-ranged interactions. The special case of vanishing 
spontaneous curvature has been investigated in Rcf. |41|. It was found that for the Canham-Hclfrich model, the 
elastic bulk and shear moduli vanish at T = 0. Small higher order curvature terms make these moduli finite, but 
thermal long-wavelength fluctuations with large amplitudes should remain. An extension of this type of analysis to 
systems with spontaneous curvature has not been attempted so far. Since the spontaneous curvature introduces a new 
length scale, many of the results of Ref. should not apply in this case. However, large fluctuations on long length 
scales can still be expected. The contribution of the fluctuations to the free-energy density should be of the form 
ksT \n-[5 / a) / , where a is the lattice constant, which is determined by the concentrations, and 5 is a molecular length 
of the amphiphile 41 4^. Therefore, the thermal contributions to the free energy of different bicontinuous cubic 
phases should be very similar, and the sequences of stable bicontinuous cubic phases as a function of concentrations 
— as shown in Fig. |^ and explained in terms of their geometrical properties in Sec. ^ — should be observed in 
experiments. In fact for the system DDAB-water-styrene a sequence of double structures was reported which seems 
to correspond to increasing topology index Js^ . For the non-cubic phases, thermal fluctuations correspond mainly to 
steric interactions between fluctuating lamellae and cylinders and translational entropy for spheres. Extending the 
model by these contributions would in fact remove the non-concavity of their free energy densities. Previous work 
suggests that the main effect for surfactant systems would be to favor micellar phases near the water-amphiphile side 
of the phase diagram and lamellar phases near the water-apex . Such a modiflcation could bring the calculated 
phase diagrams of Fig. ^ in good qualitative agreement with experimental phase diagrams measured for the whole 
Gibbs triangle, like the ones obtained in Ref. for H2O / Ciq / Ci2Ez. For a more detailed comparison, experimental 
structure determination for bicontinuous cubic phases in ternary surfactant systems is needed. The same holds true 
for ternary lipid systems where to our knowledge hardly any experimental data is known. 

With respect to long-ranged interactions, it seems reasonable to assume similar effects for all bicontinuous cubic 
phases. For attractive/repulsive forces, we therefore expect all of them to be favored/disfavored to a similar degree. 
In principle such additional contributions could stabilize the bicontinuous cubic phases with negative curvature index. 
In regard to the non-cubic phases, it was shown previously that van der Waals interactions should show a considerable 
effect only for large values of the bending rigidity, i.e. for lipids Thus we predict that the geometrical arguments 
presented in this work describe the main physics of bicontinuous cubic phases in ternary amphiphilic systems with 
spontaneous curvature, even when more complicated models are considered. 
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2.3451 
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0.4965 
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0.3173 


C(P) 


-16 
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0.6560 
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0.0466 


0.1226 


0.1636 



TABLE L Properties of the families of triply periodic surfaces of constant mean curvature, which are related to their minimal 
surface members, x ''■nd are the Euler characteristic and the surface area in the conventional unit cell (the value of Aq is 
known exactly in terms of elliptic functions for G, D, P, C(P) and I-WP |^,|8)). Fo = (^o^/27r|x|)^''^ is the topology index; the 
structures are ordered with with respect to decreasing magnitude of Fq. vq is the volume fraction of one of the two labyrinths; 
the volume fraction of the other one follows as 1 — vo- c = dv{H*) / dH*\H*^Q, where v{H*) is the volume fraction of one of the 
two labyrinths for the corresponding family of surfaces of constant mean curvature. The values for D, FWP and P are taken 
from Ref. jl^, the ones for G, F-RD and C(P) are obtained from our spline interpolation of the numerical data of Ref. [ ]l2[ . 
c' = —Ao'^/2-KX is the estimate of c as derived in the text. Note that cAq gives a very similar hierarchy as c' Aq = Fq^. 
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labyrinth 1 


surface 


labyrinth 2 
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water 


monolayer 


oil 


Q/ 


oil 


water-filled bilayer 


oil 


Qii 


water 


oil-filled bilayer 


water 



TABLE IL Distribution of oil, water and amphiphile in the different structural types of bicontinuous cubic phases. Q denotes 
single structures, Qj double structures of type I (oil-in- water) and Q// double structures of type II (water-in-oil or inverse). 
Although double structures consist of two amphiphilic monolayers, for this classification we consider each structure to consist 
of one triply periodic minimal surface separating two percolating labyrinths. 
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FIG. 1. Correspondence between (v,w)-plane and Gibbs triangle. Only a certain part of the (w, ui)-plane is mapped onto the 
Gibbs triangle. The mapping depends on the values for the amphiphile chain length q and for the spontaneous curvature cq. 
Throughout this work we use q = 1/2 (length is measured in units of amphiphile length). For co = 1/6 and Co = 1/12, the 
darkly and lightly shaded parts are mapped onto the Gibbs triangle, respectively. Water, amphiphile and oil apex of the Gibbs 
triangle are denoted by W, A and O, respectively. 

FIG. 2. Minimal surface member for the following families of triply periodic surfaces of constant mean curvature: (a) D, (b) 
C(P), (c) I-WP and (d) F-RD. Shown is one conventional unit cell. For single structures with vanishing mean curvature, these 
surfaces represent the oil-water interfaces. 

FIG. 3. Geometrical data for triply periodic surfaces of constant mean curvature: (a) scaled surface area A* and (b) scaled 
mean curvature H* as a function of volume fraction v. Each family exists only for a certain u-interval; it can be used in (b) to 
identify the family corresponding to each curve. We only show the two branches connected by the minimal surface member; 
there exist other branches within the same ii-intervals which correspond to dense arrangements of nearly spherical regions 
connected by small necks. They are assumed to have no physical relevance. The two branches used are symmetrical for G, D, 
P and C(P) since their minimal surface members are balanced. 

FIG. 4. (a) Single and (b) double structure of the G-family. The two gyroid structures are the most stable bicontinuous 
cubic phases for negative saddle splay modulus k. For the double gyroid of type I, the water forms the sheet-like region between 
the two monolayers which effectively form a bilayer in (b). For the double gyroid of type II, it fills the two channel networks 
which are separated by this bilayer. 

FIG. 5. (a) Curvature index A and (b) topology index F as a function of volume fraction v for all structures considered. 
Although the different structures can be identified according to their u-interval of existence (compare Fig. |^ , here we intend to 
demonstrate only the difference between the three structural types. For Q (solid lines), Q/ (dashed lines) and Qu structures 
(dotted lines) the curvature index vanishes and the topology index attains its maximal value for v = vq, v — 1 and « = 0, 
respectively. Qu structures and Q structures with v > vo have negative curvature indices and therefore are not stable for 
Co > 0. 

FIG. 6. Lines of vanishing frustration of various structures for r = 0. In the {v, u))-plane, the lines w{v) correspond to values 
for {w, v) where a specific phase has vanishing free-energy density. For w ^ 1 the Q and Qi structures are stable for v < vq 
and V < 1, respectively. The lines at u; = 3 and w = 4 correspond to spheres and cylinders, respectively. For cq — 1/6, only 
the part below the dotted line is mapped onto the Gibbs triangle. 

FIG. 7. Phase diagram as a function of w and r for (a) v = 0.1, (b) v — 0.2, (c) v = 0.4 and (d) v — 0.6. For co = 1/6, only 
the region below the dashed line is mapped onto the Gibbs triangle. The most stable phase with respect to r is the double 
gyroid Gj which however cannot exist for v < 0.112. 

FIG. 8. Phase behavior in the Gibbs triangle for (a) r = 0.01, cq = 1/6, (b) r = 1/15 = 0.067, co = 1/6 and (c) r = 1/15, 
Co = 1/12. For larger values of r, the gyroid structures G and Gi dominate. Lowering spontaneous curvature cq corresponds 
to raising temperature and extends the region of stability for the lamellar phase L. For simplicity, in these plots we do not 
consider the close-packing constraint for S. 
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